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ABSTRACT 

We present and discuss five candidate exoplanetary systems identified with the Kepler space- 
craft. These five systems show transits from multiple exoplanet candidates. Should these objects 
prove to be planetary in nature, then these five systems open new opportunities for the field of 
exoplanets and provide new insights into the formation and dynamical evolution of planetary 
systems. We discuss the methods used to identify multiple transiting objects from the Kepler 
photometry as well as the false-positive rejection methods that have been applied to these data. 
One system shows transits from three distinct objects while the remaining four systems show 
transits from two objects. Three systems have planet candidates that are near mean motion 
commensurabilities — two near 2:1 and one just outside 5:2. We discuss the implications that 
multitransiting systems have on the distribution of orbital inclinations in planetary systems, and 
hence their dynamical histories; as well as their likely masses and chemical compositions. A 
Monte Carlo study indicates that, with additional data, most of these systems should exhibit 
detectable transit timing variations (TTV) due to gravitational interactions — though none are 
apparent in these data. We also discuss new challenges that arise in TTV analyses due to the 
presence of more than two planets in a system. 

Subject headings: planetary systems — Stars Individual (KIC 8394721, KIC 5972334, KIC 10723750, 
KIC 7287995, KIC 7825899) — techniques: spectroscopic, photometric 
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1. Introduction 

The discovery of dozens of transiting planets 
has enabled astronomers to characterize key phys- 
ical properties of the planets, including their sizes, 
densities, atmospheric composition, thermal prop- 
erties, and the projected inclination of the orbit 



siting planet (s) (Leger et al. 2009 Queloz et al. 



with respect to the stellar spin axis (Charbon- 



neau et al. 20071. Ground-based transit searches 



have surveyed many more stars than radial ve- 
locity planet searches, allowing them to discover 
relatively rare planets, such as giant planets with 
orbital periods of less than two days. However, 
ground-based transit surveys are only efficient for 
large planets with relatively short orbital periods. 
These strong detection biases and the likely dy- 
namical instability of a system with multiple gi- 
ant planets packed close to the host star, may ex- 
plain why ground-based transit surveys have yet 
to detect a system with multiple transiting planets 
orbiting the same star. 

The Kepler mission was designed to detect 
terrestrial-size planets in the habitable zone of the 
host star, necessitating both a large sample size 
and sensitivity to a much larger range of orbital 



separations than ground-based surveys ( Borucki et 



al. 2010a). The instrument is a differential pho- 



tometer with a wide (105 square degrees) field-of- 
view (FOV) that continuously and simultaneously 
monitors the brightness of approximately 150,000 
main-sequence stars. A comprehensive discussion 
of the characteristics and on-orbit performance of 
the instrument and spacecraft is presented in |Koch| 
eraD(l20T0|. 



Its sensitivity to small planets over a wide 
range of separations gives Kepler the capability of 
discovering multiple planet systems. For closely 
packed planetary systems, nearly coplanar sys- 
tems, or systems with a very fortuitous geomet- 
ric alignment, Kepler is likely to detect transits of 
multiple planets. For systems with widely spaced 
planets or large relative inclinations, not all plan- 
ets will transit, but some may still be detectable 
based on transit timing variations (TTVs) due 
to the gravitational perturbation of one or more 



non-transiting planets (Agol et al. 2005 Holman 



and Murray 2005). In other cases, non-transiting 



planets may be detectable by follow-up observa- 
tions, such as radial velocity observations origi- 
nally intended to measure the mass of the tran- 



2009) 



Radial velocity (RV) surveys have shown that 
giant planets often reside in multiple planet sys- 
tems jWri ght et al.|200"9 ). Given the large number 
of candidate planets identified by Kepler ( Borucki 
et al. 2010a), it is expected that some fraction 



of them will be in multiple planet systems and 
a fraction of those will have multiple planets that 
transit. In addition to the ability to character- 
ize physical properites of each transiting planet, 
planetary systems with multiple transiting plan- 
ets present several advantages. For example, the 
fact that each planet formed from the same proto- 
planetary disk provides more powerful constraints 
for models of planet formation and orbital migra- 
tion. Moreover, these systems are quite powerful 
for studying the detailed orbital dynamics through 
transit timing variations (or lack thereof) . In some 
cases, the planet masses may be determined with- 
out measurements of stellar RV variations. Sys- 
tems with more than one transiting planet hold 
unique power in this regard as the period, orbital 
phase, and approximate size of the various planets 
are known. This information can help significantly 
in finding a unique solution to an otherwise chal- 



lenging and degenerate inversion problem (Ford 



and Holman||2007| [Nesvorny and Morbideni|[2008 



Meschiari and Laughlin||2010 1 

In cases where RV measurements are able 
to measure the planet masses independently of 
TTVs, the two techniques can be combined to 
measure the mass and size of the host star with- 
out relying on stellar models (see Agol et al.|2005 
Holman and Murray 2005). In cases where RV 



observations are not practical (e.g., hot stars, fast 
rotators) or would require prohibitive observing 
time (i.e., faint stars), the detection of TTVs can 
be used to confirm that transit candidates orbit 
the same star — as opposed to being two objects 
transiting two stars blended within a single point- 
spread-function (PSF) — and to determine if the 
companions are of planetary mass. Here multi- 
ple transiting systems are particularly powerful 
as the period, phase, and size of additional plan- 
ets can be determined from the light curve. This 
additional information can also help simplify the 
inverse problem (Steffen and Agol||2007[ ) . 



We present five planetary candidate systems in 
which the transits of multiple objects can be seen 
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in the first quarter of photometric data (a 33.5-day 
data segment from May 13 to June 15 UT, 2009) 
from the Kepler spacecraft. While not confirmed 
planet discoveries, these systems have passed sev- 
eral important tests that eliminate false-positive 
signals. If all were ultimately shown to be plan- 
ets, then these systems would contain four plan- 
ets with radii smaller than three Earth radii (the 
smallest being two Earth radii), at least two pairs 
of planets in or very near a low-order mean-motion 
resonance (MMR), and one system with at least 
three distinct transiting planets. 

For simplicity, we will refer to these objects as 
"planets" throughout this paper, recognizing that 
their confirmation as such is yet incomplete and 
that some of these transit signals may be due to 
other astrophysics. The stellar references that we 
will use throughout this paper are Kepler Objects 
of Interest (KOI) 152, 191, 209, 877, and 896 with 
the transiting planets denoted by ".01" , ".02" , etc. 
beginning in the order that they were identified 
with the transit detection software from the Ke- 
pler pipeline. Thus, the planet number designa- 
tion does not necessarily reflect the order of the 
planets within each system. We do not use letter 
designations, which by convention are reserved for 
confirmed planets. 

This paper will proceed as follows. First, we 
give the known properties of the host stars (Sj2|. 
In ^ we discuss the photometric reduction and 
the algorithm used to identify the multiple can- 
didates within each system. We also outline the 
tests we have conducted to eliminate false-positive 
systems. We present estimates of the orbital and 
physical properties of these objects should they 
prove to be planets (Q. In ^j5]we discuss the pos- 
sible future detection of transit timing variations 
based upon a Monte Carlo simulation of these can- 
didate systems. Finally, we discuss the implica- 
tions of these results in J|6j 

2. Stellar properties 

For each of the five stars, we obtained high res- 
olution echelle spectroscopy from the McDonald 
Observatory 2.7m telescope Tull Coude spectrom- 
eter with resolving power R=60,000. We also ob- 
tained one spectrum of KOI 155 with the FIES 
spectrograph on the Nordic Optical Telescope and 
one spectrum of KOI 191, KOI 209, KOI 877, and 



KOI 896 with the Kitt Peak National Observatory 
4-m telescope. These spectra were obtained for the 
purpose of constraining effective temperature T^g, 
surface gravity logg, projected rotational velocity 
for the star vsini, and metallicity [M/H]. The 
McDonald spectra were reduced and extracted us- 
ing the IRAF echelle package 

In all cases, the spectroscopic analysis was done 
by matching the observed spectra to a library of 
synthetic spectra. The synthetic spectra cover the 
wavelength region 5050-5360 A (centered roughly 
on the Mg b lines). The grid has coarseness of 
250 K in T^g, 1-4 km/s in u^ot — wsini, and 0.5 
dex in logy and [M/H], implying uncertaintes of 
half those values. For the host stars with low SNR 
spectra, we performed a diagnostic using the J— K 
color to verify that their compositions are con- 
sistent with solar metallicity, but the values we 
report are simply from template-matching with 
[M/H] fixed at 0. Table [ijhsts the resulting stel- 
lar parameters for all five stars as well as the in- 
struments used in the observations. All five stars 
apparently reside near or on the main sequence. 

3. Kepler data and photometric analysis 
3.1. Transit identification 

Each of these systems was found using the 
Transiting Planet Search Pipeline (TPS) which 
identifies significant transit-like features, or Thresh- 
old Crossing Events (TCE), in the Kepler light 
curves ( Jenkins et aL|[2010 ). Data showing TCEs 
are then passed to the Data Validation (DV) 
pipeline ( |Wu et al.\20Tol . The purpose of the DV 
pipeline is twofold: 1) to fit a transiting planet 
model to the data, remove it from the light curve, 
and to return the result to TPS in an effort to 
find additional transit features and 2) to com- 
plete a suite of statistical tests that are applied 
to the data after all TCEs are identified in an 
effort to assess the likelihood of false-positives. 
The binary discrimination statistics and the mo- 
tion detection statistic, in particular, speak to 
the likelihood of astrophysical false positives such 
as grazing eclipses and diluted eclipsing binaries. 
These statistics are described below. 



^IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Uni- 
versities for Research in Astronomy (AURA) under coop- 
erative agreement with the National Science Foundation. 
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Table 1 

Stellar properties and locations for the five candidate systems. 



KOI 


KIC-ID 


ra 


DEC 


Kcpmag 






[M/H] 


V sin i 










(2000) 


(2000) 




(K) 






(km/s) 


(km/s) 




152 


8394721 


20 02 04.1 


44 22 53.7 


13.9 


6500 


4.5 


(fixed) 


14 


-22.32 




191 


5972334 


19 41 08.9 


41 13 19.1 


15.0 


5500 


4.5 


(fixed) 





-62.98 


1" 


209 


10723750 


19 15 10.3 


48 02 24.8 


14.2 


6100 


4.1 


-0.05 


7.8 


-12.78 


4" 


877 


7287995 


19 34 32.9 


42 49 29.9 


15.0 


4500 


4.0 


(fixed) 





0.531 


2t,,c 


896 


7825899 


19 32 14.7 


43 34 52.9 


15.3 


5000 


4.0 


(fixed) 


1 


-21.28 


1^ 



Note. — Each system was analyzed by matching the observed spectra to the CfA hbrary. 
Telescopes used for these observations: 
a, Nordic Optical Telescope, FIES 
6, McDonald Observatory 2.7m, Tull Coudo 
c, Keck-HIRES 



After pipeline data processing and the photom- 
etry extraction, the time series is detrended with 
a running 1-day median. All observations that oc- 
cur during transit are not included in the evalua- 
tion of the median. The transit lightcurve is mod- 



eled using the analytic expressions of Mandel and 



Agol ( 2002 1 using the non-linear limb darkening 



parameters that Claret ( 2000 ) derived for the Ke- 
pler bandpass. A first estimate for Af* and i?* is 
obtained by comparing the derived stellar T^g and 



log g values, obtained using the method of ( Valenti 



and Piskunov 1996), to a set of CESAM (Morel 



1997) stellar evolution models computed in steps 
of 0.1 Mq for solar composition. We note that in 
some cases, particularly KOI 877 which shows sig- 
nificant spot modulation, there is some residual 
noise power at time scales longer than the tran- 
sit duration. While this noise causes small base- 
line fluctuations in the out-of-transit portion of 
the light curves, it does not significantly affect es- 
timates of the transit model parameters (e.g., the 
transit depth). 

With and i?* fixed to their initial values, 
a transit fit is then computed to determine the 
orbital inclination, planetary radius, and depth 
of the occultation (passing behind the star) as- 
suming a circular orbit. The best fitting model 
is found using a Levenberg-Marquardt minimiza- 
tion algorithm (Press et al. 1992). This model is 



then removed from the lightcurve and the residu- 
als are used to characterize the next transit candi- 
date identified by TPS. The light curves and tran- 
sit models for these five systems are shown in Fig- 
ures[TJ[2]|3j|4j and[5] The significance of individual 



transit events for each planet and the significance 
of the folded transit event (in terms of canonical 
"sigmas") are given in the first two columns of 
Table [21 Note that these statistics come from the 
pipeline processing and are not calculated for sys- 
tems where there is only one transit of the most 
significant transit event (KOIs 152 and 209); the 
additional candidates in these two systems were 
found manually. 

3.2. Candidate vetting from Kepler data 

The depths of planetary transits should be 
consistent from orbit to orbit as well as evenly 
spaced in time. Eclipsing stellar binaries, on 
the other hand, generally have primary and sec- 
ondary eclipses with different depths that are fre- 
quently spaced asymmetrically in time. The bi- 
nary discrimination test includes two metrics: the 
Odd/Even statistic and the Epoch statistic. The 
Odd/Even statistic is a comparison between the 
depth of the phase-folded, odd-numbered transits 
and the depth of the phase-folded, even-numbered 
transits. The Epoch statistic compares the tim- 
ing of the odd and even-numbered transits. Both 
statistics are constructed as distributions, and 
the significance (reported in Table [2]) of the statis- 
tic is obtained by evaluating the cumulative dis- 
tribution function for the appropriate number of 
statistical tests. The significance is the probabil- 
ity that the statistic is consistent with the binary 
interpretation. 

The motion detection statistic identifies objects 
with flux-weighted centroids that are highly corre- 
lated with a transit signature derived from the flux 
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Fig. 1. — Unbinned light curve and transit models 
for the candidates KOI 152.01 (top), 152.02 (mid- 
dle), and 152.03 (bottom). The lower curves in 
each panel show the transit while the vertically- 
offset, upper curves show the data at phase = 
0.5. The "+" and "*" symbols are for even and 
odd transits. The top-most panel shows the raw 
lightcurve with the transits for the three candi- 
dates marked with vertical lines. 
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Fig. 2. — Light curve and transit models for the 
candidates KOI 191.01 (top) and 101.02 (bottom). 
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Fig. 3. — Light curve and transit models for the 
candidates KOI 209.01 (top) and 209.02 (bottom). 



Fig. 4. — Light curve and transit models for the 
candidates KOI 877.01 (top) and 877.02 (bottom). 
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Fig. 5. — Light curve and transit models for the 
candidates KOI 896.01 (top) and 896.02 (bottom). 



timeseries. The statistic is a variable with two 
degrees of freedom (row and column). The prob- 
ability of producing a statistic of equal or lesser 
value is reported as "Motion" in Table [2j It is 
computed by evaluating the cumulative distri- 
bution function at the value of the detection statis- 
tic given two degrees of freedom. The complement 
of this value is reported so that values near unity 
represent a small likelihood of a correlation. Cen- 
troids that are highly correlated with the transit 
signature are indicative of a crowded photometric 
aperture — a warning that the transit could be due 
to a nearby eclipsing binary diluting the flux of 
the target star. Such cases are denoted by a mo- 
tion significance near zero. A full description of 



DV statistics is given in Wu et al. (2010). 



Outside of the pipeline, all TCEs are fitted with 
a planet transit model as described in ^Batalha ^ 



al. (2010). Those yielding an estimated planet ra- 
dius less than 2 Rj are assigned a Kepler Ob- 
ject of Interest (KOI) number. The modeling re- 
turns an independent test of the Odd/Even statis- 
tic (Depth-sig in Table [2]), expressed in units of 
the standard deviation. The modeling also tests 
for the presence of secondary eclipses (or occulta- 
tions) at phase = 0.5 and reports this as Eclipse- 
sig in Table [2j also in units of the standard devia- 
tion. We note that the transit modeling "Depth- 
sig" statistic and the DV "Odd/Even" statistic are 
not identical as they arise from different analyses 
(especially when the odd and even transit depths 
are not significantly different) but will be approx- 
imately equivalent when a significant detection is 
seen. 

Figure [6] shows the normalized relative flux vs. 
centroid position (referred to as a "rain plot") as 
described in |Batalha et al. (2010) for each target 
star. Here, the relative flux is plotted against the 
relative centroid position along rows and columns. 
A centroid shift that is highly correlated with the 
transit signature would appear as a diagonal de- 
viation in the plot (see bottom-right panel of Fig- 
ure[6]) whereas uncorrelated centroids "rain down" 
vertically under the out-of-transit points. No sig- 
nificant correlations can be readily identified in 
the rain plots for any of the candidates presented 
here except KOI 191. This is reflected in the Mo- 
tion statistics reported in Table [2] KOI 191 has 
a correlation with a 98% signiflcance (taking the 
complement of the value reported in Table |2| — 
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Table 2 

Planetary candidate transit detection and validation statistics. 



KOI 


Transit-Sig 
(Single) 


Transit-Sig 
(Folded) 


Odd/Even 


Epoch 


Motion 


Depth-Sig 


Eclipse-Sig 


152.03 












0.6 


-0.5 


02 












1.3 


-1.0 


01 
191.02 


4.2 


13.0 


0.71 


0.95 


0.56 


0.0 


-0.5 


01 


104 


211 


0.43 


0.97 


0.02 


1.8 


-2.6 


209.02 












0.3 


-0.7 


01 
















877.01 


4.9 


13.0 


0.14 


0.97 


0.93 


1.7 


2.5 


02 


5.0 


10.4 


0.83 


0.97 


0.94 


0.5 


0.5 


896.02 


5.3 


13.3 


0.11 


0.96 


0.61 


0.0 


-4.0 


01 


7.9 


15.7 


0.16 




0.45 


1.9 


1.5 



comparable to a 2.3-sigma detection. A correla- 
tion does not necessarily rule out the planetary 
interpretation, rather it should be interpreted as a 
warning that the photometric aperture is crowded. 
Additional analysis or observational follow-up is 
required to determine the location and magnitude 
of each star in the vicinity and, ultimately, the 
origin of the transit-like features. 

None of the DV binary discrimination statistics 
reported in Table [2] are significant at the 3-sigma 
level or higher (corresponding to a significance of 
0.997 or larger). However, DV statistics for KOI 
152 and 209 are not available due to the fact that 
the dominant transit feature in the Quarter 1 light 
curve appears only once. A planet transit model 
cannot be fitted to a light curve with a single tran- 
sit with DV. Consequently, it does not get filtered 
and passed back to TPS for the detection of the 
shallower, shorter-period transits. The shorter- 
period events are subjected to light curve mod- 
eling and the associated statistics are reported. 

3.3. Candidate vetting from ground-based 
data 

In order to identify neighboring stars that are 
located within a few arcsec of the target star, we 
obtained optical images of all five stars using the 
photometries CCD guide camera on the HIRES 
spectrometer on the Keck 1 telescope (shown in 
Figure^. The images were obtained with no filter 



over the CCD that has high sensitivity from 400- 
800 nm, mimicking the CCD on the Kepler cam- 
era. The images were taken in seeing of 0.7 - 0.9 
arcsec and clear skies. Each image is 43 x 57 arc- 
sec, with 0.30 arcsec/pixcl. Such neighboring stars 
bring two concerns, diluting the transit depths and 
possibly being eclipsing binaries that cause the 
transit signal. In the latter case, one or more of 
the candidate planets may be false positives. 

For KOI 152, the Keck image shows two neigh- 
boring stars located south and east of the KOI 
(which is the brightest star in the field). The 
brighter of the two neighbors is 5.4 arcsec SE and 
3.8 mag fainter. The fainter neighboring star is 
4.5 arcsec to the southeast and 6.6 mag fainter. 
The Keck image of KOI 209 reveals no neighbor- 
ing stars down to 20th mag. The Keck image of 
KOI 896 reveals one neighboring star 7.3 arcsec to 
the southeast that is 2.9 mag fainter, and another 
neighboring star 8.1 arcsec to the north that is 1.5 
mag fainter. The wings of both neighboring stars 
encroach into the Kepler photometric aperture. 
The Keck image of KOI 877 shows no neighboring 
stars down to 20th mag. We have not examined 
these neighboring stars to determine if they are 
eclipsing binaries. However, the centroid statis- 
tics given above indicate that, should they prove 
eclipsing binaries, they are unlikely to be the cause 
of the observed transit signature. 

For KOI 191, the Keck image reveals a neigh- 
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Fig. 6. — Relative flux vs. centroid position for 
KOIs 152, 191, 209, 877, and 896 beginning in 
the top left corner. The plot in the bottom right 
is from a background eclipsing binary and is pre- 
sented to illustrate the difference in the data be- 
tween good planetary candidates and background 
eclipsing binaries. The blue squares correspond 
to rows while the red circles are columns on the 
CCD. 



boring star located 1.5 arcsec east of the main 
star and 2.6 mag fainter. As stated above, KOI 
191 exhibits some correlation in the rain plots 
which indicates a crowded field. The light curve 
from KOI 191 also shows an additional, periodic 
transit feature. The ephemeris of this feature is 
Tc - 2454900571)) = 65.6589 + Ex 0.7086 Days. 
Its V-shaped transit shape (0.2 mmag depth) and 
its two-hour duration indicate that a faint eclips- 
ing binary is also blended with the target. 

3.4. Blender analysis 

False positive scenarios were investigated for 
the five systems by exploring the possibility that 
the Kepler photometry is the result of contamina- 
tion of the light of the candidate by an eclipsing bi- 
nary along the same line of sight, a "blend" . Given 
that the centroid motion statistics discussed above 
rule out a large fraction of the background con- 
taminants, we focused here on hierarchical triple 
systems in which the candidate and the binary are 
at the same distance. Angular separations in these 
cases would usually be too small to generate sig- 
nificant centroid motion. 

We modeled the Kepler photometry of each 
candidate assuming it is the result of the bright- 
ness variations of an eclipsing binary being atten- 
uated by the (typically) brighter candidate star. 
For KOIs with two or more signals in the light 
curve, we modeled the light curves at each pe- 
riod separately and accounted for possible blends 
at the other period(s) by incorporating extra di- 
lution consistent with those other stars. In these 
cases, whether the blended eclipsing systems at 
each period are related or not (i.e., in a hierar- 
chical quadruple system, for candidates with two 
signals) is immaterial for the purposes of model- 
ing the light curves. The objects composing the 
binary are referred to as the "secondary" and "ter- 
tiary", and the candidate is the "primary". The 
procedure closely follows that described by [Torres] 
et al. (2004) and consists of calculating synthetic 



light curves that result from the three objects for 
a wide range of eclipsing binary parameters, and 
comparing those light curves against the Kepler 
photometry for the candidate, in a sense. We 
regard as acceptable any blend scenario that re- 
sults in a synthetic light curve giving a for the 
fit that is not significantly different (at the 3-ct 
confidence level) from a planet model fit. 
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Fig. 7. — Optical images of all five Kepler stars 
taken in seeing of 0.7-0.9 arcsec and clear skies 
by the guider camera on the Keck-HIRES instru- 
ment. The images reveal neighboring stars that 
may affect the interpretation of the transits, no- 
tably their depths and the existence of eclipsing 
binaries. In order from the top left the images are 
of KOIs 152, 191, 209, 877, and 896. 



The brightness variations of the binary are gen- 
erated using detailed calculations including limb 
darkening, gravity brightening, reflection, and 
proximity effects. The properties of the candi- 
date are tightly constrained by the spectroscopic 
parameters in Table [l] and were held fixed. The 
parameters of the binary components were taken 



from model isochrones by Girardi et al. (20001, 



parametrized in terms of their mass. The sec- 
ondary and tertiary masses were allowed to vary 
over wide ranges (0.1-1.4 Mq) in order to fit the 
light curve, and the inclination angle was also a 
free parameter. By comparing the quality of the 
light-curve fits over these ranges, we are able to 
constrain the properties of the secondaries and 
tertiaries that provide acceptable fits, and we can 
reject other blends. We consider hierarchical con- 
figurations of two types: ones in which the tertiary 
is a star, and ones in which it is a planet contribut- 
ing no light. In the latter case the size of the planet 
is a free parameter, which we varied between 0.1 
and 2.0 i?jup- We account for the additional stars 
identified by high resolution imaging by including 
the proper amount of extra light in our models. 

In all five cases we find that configurations with 
stellar tertiaries are inconsistent with the Kepler 
light curves, for any size secondary. Thus, hierar- 
chical triples involving a stellar binary are ruled 
out. When the tertiary is allowed to be a planet, 
we find that there is a range of possible solutions 
that yield acceptable matches to the Kepler pho- 
tometry, often times as good as obtained from 
a single star and planet model. However, many 
of those solutions can also be excluded on other 
grounds as described below. 

For KOI 152, which is a mid F dwarf, the only 
configurations consistent with the Kepler light 
curves involve secondary stars that are at least as 
massive as the primary — and therefore almost as 
bright or brighter. In all cases the planetary com- 
panions are roughly 0.5 i?jup in size. Such bright 
secondaries are unlikely as they would have been 
seen spectroscopically, unless the velocity differ- 
ence with the target star happens to be very small 
so that the lines are indistinguishable. This can- 
not be completely ruled out if the two stars are 
in a wide orbit around each other, but it would 
require a special set of circumstances given the 
constraints on the centroid motion. Thus, except 
for this particular case of a bright secondary, our 
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analysis supports the planetary interpretation. 

For KOI 209, a mid to late F dwarf, we again 
find that the only blend configurations that fit the 
Kepler light curve involve a secondary that is as 
bright as the primary or brighter, and a tertiary 
about the size of Jupiter. A planetary interpreta- 
tion is again favored. 

Our blend simulations for KOI 896, an early K 
star, again indicate that the secondaries required 
for a good fit to the light curve are similar in 
brightness to the primary, and the tertiaries have 
sizes between 0.4 and 0.5 -Rjup- Solutions with 
smaller secondary stars are visibly worse. There- 
fore, except for the somewhat artificial case of twin 
stars, the Kepler photometry favors a planetary 
interpretation. 

For KOI 191.02 the light-curve fits allow blend 
scenarios in which the secondaries are smaller than 
the late G-type primary, down to a spectral type 
of late K. The tertiaries tend to have radii around 
0.3 i?jup. In the case of KOI 191.01 we find 
that the light curve prefers secondaries that are 
as bright as or brighter than the primary and ter- 
tiaries of '--^l.S i?jup- Both of these results take 
account of the fact that there is extra light from 
the 1.5 arcsec companion described earlier, which 
we have considered here to be an unrelated back- 
ground eclipsing binary. 

KOI 877 is a late K dwarf. In addition to the 
blend solutions with bright secondary stars, we 
find acceptable fits to the light curves with sec- 
ondaries up to two magnitudes fainter than the 
primary (spectral type M2-M3), which might not 
be noticed spectroscopically. The size of the ter- 
tiaries in these cases is about 0.4 i?jup. To rule out 
such configurations, additional observational con- 
straints are needed, such as accurate multi-band 
photometry out of transit to check for color incon- 
sistencies (e.g., O'Donovan et al.||2006 ). 

We note that two eclipsing binaries diluting 
each other, perhaps in a hierarchical quadruple 
system (e.g., Pilecki and Szczygiel|2007 ) is among 
the most likely, still-viable false positive scenar- 
ios that involves no planets. To our knowledge, 
there are no well-established cases of stellar sys- 
tems in which more than one body passes in front 
of the primary. The All-Sky Automated Survey 



ries" : physically bound quadruples or blended in- 
dependent binaries in which each pair of binaries 
eclipses. Such systems are viable false positive sce- 
narios since they produce no astrometric centroid 
motion, can be dynamically stable, and are as- 
trophysically reasonable. (Two small stars orbit- 
ing the same primary can be dismissed — except 
perhaps for KOI 191 — because at the periods ob- 
served here, such systems are not stable even on 
very short timescales.) The probability that all 
of the transiting planet candidates in a multiple 
system are produced by eclipsing binary false pos- 
itives is less than ~ 10^'^ since roughly this frac- 
tion of the Kepler stars are eclipsing binaries ( Prsa 



et al. 2010) and the two must be blended into 



a single PSF. Similar arguments, using the ap- 
proximate number of Kepler planetary candidates 



(Borucki et al. 2010b), show that the scenario of 



(Pilecki and Szczygiel 2007) identified seven sys- 



tems that were termed "double eclipsing bina- 



two independent stars with independently orbiting 
planets must be smaller still. 

In summary, the above blend analyses suggest 
that KOI 152, 209 and 896 are very likely true 
planets orbiting a common star. KOI 191 and 877, 
on the other hand, could be blends of two planets 
orbiting two separate stars. In none of these sys- 
tems was a configuration of transiting stars pre- 
ferred over transiting planet-sized objects. 

4. Planetary and system properties 
4.1. Orbital properties 

Of the five candidate planetary systems, there 
are four pairs of planets with a well-characterized 
orbital period ratio. Two pairs are near the 2:1 
MMR (KOI 152.02/01 and KOI 877.01/02), one 
pair is slighly outside the 5:2 commensurability 
(KOI 896.02/01), and one pair is hierarchical with 
a period ratio exceeding 6:1 (KOI 191.02/01). The 
proximity of KOI 152 and 877 to the 2:1 MMR 
hints that these systems are likely to have large 
TTVs. However, the timescale for large TTVs in 
resonant systems can be quite long. For a sys- 
tem of Neptune-mass planets librating about a 2:1 
MMR, we would not expect to detect TTVs based 
on the QO and Ql data presented. Moreover, for 
librating systems with smaller masses or that are 
very close to exact resonance, the time needed to 
distinguish a TTV signal from a constant period 
lengthens, also requiring additional data. 

Inferring the relative frequency of planets with 
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various orbital spacings will require a population 
analysis that corrects for the geometric transit 
probabilities. Of course, there may be additional 
non-transiting planets, or small transiting plan- 
ets, between KOI 191.02 and 191.01 or in any of 
the candidate systems presented. In some cases, 
TTVs or radial velocity follow-up could identify 
such planets. In other cases, the non-transiting 
planets will remain undetected, further complicat- 
ing a population analysis. The sample of multi- 
ple planet candidate systems presented here is too 
small — and not necessarily unbiased — for such an 
analysis. 

For two of the transit candidates presented, 
only one transit has been observed. There is a 
small chance that a second transit occured during 
a data gap. However, this is a priori unlikely and 
the extended transit durations also support large 
orbital periods. The periods listed in Table [4T] are 
lower limits based upon the non-observation of a 
second transit in the data. Also included, however, 
are estimates based upon the transit duration as- 
suming a circular orbit and a central transit of the 
planet. With these latter estimates, it is tempt- 
ing to identify the KOI 209 system as lying near a 
3:1 MMR and the KOI 152 system as being near 
a 4:2:1 MMR. However, we caution that the un- 
certainty in the orbital period estimated from a 
single transit is far too large to have confidence 
that these systems are near MMR. Nonetheless, 



Yee and Gaudi ( 2008 1 show that orbital period es- 



timates based upon single transit durations need 
not be as conservative as the estimates that we 
state. 

For each neighboring pair of planet candidates, 
we measure the ratio ^ = {D„J D out) {Pout/ PinY^^ , 
where D is the transit duration and P is the or- 
bital period. In each case the ratio is near unity 
(see Table [s]), as expected for a pair of objects 
on circular and coplanar orbits around a com- 
mon star. We compare these ratios to the results 
of Monte Carlo simulations of the distribution of 
^ (denoted ^mc) for pairs of planets on circular 
and coplanar orbits with orbital periods similar 
to those observed. We assume random viewing 
angles, subject to the constraint that both planets 
transit. The 16th, 50th, and 84th percentile of 
these distributions are included in Table [S] In no 
case do we find evidence for a large eccentricity or 
for a blend of multiple stars each with one tran- 



siting object. The largest deviation from unity is 
that for KOI 191, but in this case the observed 
and expected values are very close, due to the 
large ratio of semimajor axes. In the cases of 209 
and 896, the ratio is slightly less than expected for 
coplanar, circular orbits. However, this could be 
easily reconciled if one planet in each system were 
to have a modest eccentricity (~ 0.05 for KOI 209, 
- 0.1 for KOI 



4.2. System coplanar ity 

One important question surrounding these 
multi-candidate systems is the geometric prob- 
ability that Kepler would see both planets tran- 
siting. Following Ragozzine and Holman (20101, 
and based on the method of IBorucki and Sum-I 



mers 



(19841, we calculate this probability by con- 



sidering the area of the region on the celestial 
sphere, centered on the star, that is aligned to 
see both planets transit (see also Beatty and Sea- 



ger 



2010 Gillon et al. 20101. We also note here 



the analytical approximation given by | Ragozzine] 
and Holman (2010) for the probability of observ- 



ing both planets transit as a function of the true 
mutual inclination between the planets, (j). The 
result is different in low and high mutual inclina- 
tion regimes, with the critical angle between the 
regimes as (pcin — R*/ai where i?* is the radius 
of the star and oi is the orbital distance of the 
innermost, transiting planet. In the low mutual 
inclination regime, the probability of seeing both 
planets transit is i?*/a2 where 02 is the orbital 
distance of the outer planet. Thus, the probabil- 
ity of seeing both planets transit is equal to the 
probability that the more distant planet transits. 

In the high mutual inclination regime, this is no 
longer true, and only observations along the line 
of nodes of the orbital planes will see both plan- 



ets transit (Koch 1995 Holman and Murray 2005 1. 



In this regime, the probability of observing a tran- 
sit is i?^/(aia2 sin0), significantly lower than the 
probability in the low inclination regime. With 
the three candidates seen in the KOI 152 system, 
the probability is more complicated as another 
mutual inclination angle and mutual nodal angle 
are required to specify the system. Unlike in the 
two-planet case, it is easy to construct high mu- 
tual inclination systems where no observer would 
see three-planets transit. If both mutual inclina- 
tions are low (i.e., below 1.5°), then the proba- 
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Table 3 

Orbital periods and transit epochs for the candidate planets. 



Candidate 


Transit Epoch 
BJD -2454900 


Period 
(Days) 


Period Ratio 
(vs. inner) 


Transit Duration 
(Days) 


(obs.) 


Kmc 
(predicted) 


152.03 
02 
01 


69.622 ± 0.0053 
66.630 ± 0.0079 
91.747 ±0.0026 


13.478 ±0.0098 
27.406 ±0.0150 
> 27 (51.9) 


2.03 
(3.85) 


0.2071 ±0.0022 
0.2823 ± 0.0060 
0.3432 ± 0.0013 


0.9291 
1.0188 




191.02 
01 


65.50 ±0.16 
65.3847±4 X 10"'* 


2.420 ± 0.0006 
15.359 ±0.0004 


6.347 


0.0948 ±0.0016 
0.1494 ±0.0002 


1.1751 


1 1C+0.60 
^■^^-0.13 


209.02 
01 


78.822 ± 0.0046 
68.635 ± 0.0036 


18.801 ±0.0087 
> 29 (49.3) 


(2.62) 


0.2884 ±0.0018 
0.4252 ± 0.0007 


0.9429 


-I 19+0.68 
^■^^-0.11 


877.01 
02 


103.952 ± 0.0028 
114.227 ±0.0051 


5.952 ±0.0024 
12.039 ±0.0077 


2.023 


0.0962 ±0.0012 
0.1192 ±0.0021 


1.0204 




896.02 
01 


107.051 ±0.0028 
108.568 ± 0.0024 


6.311 ±0.0024 
16.242 ±0.0075 


2.574 


0.1278 ±0.0016 
0.1916 ±0.0017 


0.9144 


1 11+0.55 
^■^^-0.10 



Note. — The periods and period ratios hsted in parentheses are estimates based upon the duration of a single 
transit. The values hsted for the predicted estimates ^mc correspond to the median with the errors corresponding 
to the 16th and 84th percentiles. 



bility is P ~ R*/az — 0.017, where 03 is the or- 
bital distance of the third planet. Introducing a 
larger mutual inclination between any pair of plan- 
ets can significantly reduce this probability. Even 
if one pair of planets has a mutual inclination of 
only 10°, the probability of seeing all three transit 
drops to ~ 0.0025. Based upon these probabilities, 
if these three objects are confirmed to be multiple 
planetary systems, then they are very likely copla- 
nar to within a few degrees. 

Next, we consider the expected number of sim- 
ilar systems for which the outer planet does not 
transit. This requires a calculation of the proba- 
bility of seeing the outer planet transit given that 
the inner planet is known to transit. This can also 



be answered with the model of Ragozzine and Hol- 



man 



(2010) and also depends strongly on the mu- 
tual inclination. Instead of providing analytical 
estimates, Figure |8] shows a Monte Carlo numeri- 
cal calculation of the fraction of random observers 
that see both planets transit out of those observers 
who see the inner planet transit. All of the KOI 
systems are shown in this figure, including two 
curves for KOIs 152.02/152.01 and 152.03/152.02. 
Even in the coplanar case, for each observed 




5 10 15 

Mutual Inclinotion (degrees) 



Fig. 8. — Geometric transit probabilities, assum- 
ing that the inner planet transits, as a function of 
mutual inclination for multi-transiting systems us- 
ing the method described in the text and discussed 
in Ragozzine and Holman (2010). The 



curves 



shown are for the 4 double-transiting systems, the 
KOI 152.02/01 pair, and the KOI 152.03/02 pair. 
These calculations assume solar-like stars and ne- 
glect the size of the planet. 
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multi-transiting system there are a-ijax multiple 
systems where only the inner planet transits. Hi- 
erarchical systems like KOI 191, must have at least 
3.5 times as many counterparts where only the 
2-day period planet is seen in transit. If any of 
these systems have large mutual inclinations, the 
number of implied similar systems increases con- 
siderably. The short orbital periods and the likely 
near-coplanar state of these systems has implica- 
tions for their formation. If these planets formed 
beyond the snow line, some mechanism must be 
invoked to bring them planets to their current lo- 
cations. We note that none of these systems are 
candidates for formation by Kozai cycles with tidal 
friction due to perturbations by a distant stellar 



brown dwarf (e.g., CoRoT-3b, 21.7Mj, [Peleuil et 



companion (Fabrycky and Tremaine 2007), since 



the presence of the other planet would shut off 
Kozai effects. 

The two major classes of remaining theories for 
moving these planets in are planet-planet scat- 



tering (Rasio and Ford 19961 and disk migration 
dGoldreich and Tremaine||1980| iLin et al.||1996|. 



Planet scattering tends to excite orbital eccentric- 
ity and inclination and has difficulty migrating 
planets into short period orbits. Disk migration 
is able to migrate multiple planets into short pe- 
riod orbits and tends to damp inclination. The 
near-resonant ratios of the observed systems fa- 
vors the disk migration hypothesis. However, con- 
tinued migration after resonant trapping excites 



eccentricities and inclinations ( Lee and Peale 2002 
Lee and Thommes|[2009l ) . 



4.3. Physical properties 

We consider the physical properties of plan- 
ets (e.g., mass and composition) that can be de- 
rived from the Kepler photometry. The primary 
motivation for this effort is to identify plausible 
masses for these planets in order to conduct a 
Monte Carlo study of TTV signals (Section [5|. 
For a given radius we estimate a range of masses 
that depends heavily on the possible bulk compo- 
sitions. For planets with radii larger than Saturn, 
the planet mass is largely indeterminate because 
of the transition in the mass-radius relation from 
a Coulomb to an electron degeneracy dominated 
equation of state. The mass-radius curve turns 
over, so an object with a IRj radius could be 
anything between a sub-Saturn mass planet (e.g., 
HAT-P-12b, 0.2Mj, [Hartman eFal] pOOg])) to a 



al. (20081). For a planet radius up to that of 
Neptune, the planet mass is constrained better be- 
tween low-density objects rich in gas and volatiles 
and the rocky, iron-rich and high-density super- 
Earths. Interpreting the bulk composition of such 
planets is more difficult due to the degeneracies 
that arise with materials having different equa- 
tions of state. The measured planetary radius and 
the expected mass ranges of the candidate planets 
is shown in Table ID 

Under these limitations, we estimate a mass 
range for each object using theoretical models, 
that are consistent with this level of observational 



uncertainty (jValencia et al. 


2006 2007 


Fortney 


et al. 2007 


Seager et al.|2007| 


Grasset et al.|2009|. 



These models of planetary interiors cover a wide 
range of physical constitutions — from pure hydro- 
gen to pure iron planets. Obviously, there would 
be extremes that could not arise in nature. One 
can constrain masses based on pure iron and pure 
water super-Earths, arguing that planet formation 
scenarios would not allow for such pure constitu- 
tions ( [Valencia et al.|2007| [Marcus et al.|2010a|b |. 

In primary planet formation of any flavor, giant 
impacts and late water delivery are the only plau- 
sible way to "purify" an initially mixed-materials 
formation in a protoplanetary disk. For example, 
the iron-enhanced bulk composition of Mercury is 
explained by an early head-on impact with a sim- 



ilar body. Marcus et al. (2010a I find that a mass- 



dependent limit on final mean density (hence, ra- 
dius) should exist for super-Earth planets more 
massive than 1 Me, which is significantly less 
dense than pure iron. On the low-density bound 



(high radius), Marcus et al. (2010b) show that 
more than about 75% by mass enrichment in pure 
water is not possible, but here the upper envelope 
is not easily constrained due to the possible ad- 
dition of a H/He envelope and/or extended atmo- 
sphere for a hot planet ( Adams et al.|2008 Rogers 
and Seager][2010[ ). 

Starting with the small-size objects, KOI 
191.02, with Rp = 2.{)Re, may well be a super- 
Earth. It is half the size of the ice giant Uranus 
{ARe) and 30% smaller than GJ1214b {2.7 Re). 
However, the mass range of 5 — ISM^ spans the 
range between a water-rich world and an iron-rich 



remnant of a giant impact collision ( Marcus et al 



2010a). We use a mass of Mp = IOMe for KOI 
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Table 4 

Planetary radii, fractional error, and likely range of masses. 



Candidate 


Planet Radius 


Fractional Radius 


Mass Range 




(Rj) 


(Re) 


Error 




152.03 


0.30 


3.36 


7.8% 


9 - 30 Mb 


02 


0.31 


3.47 


6.7% 


9-30 Me 


01 


0.58 


6.50 


1.2% 


20 - 100 Me 


191.02 


0.18 


2.04 


6.6% 


5-18 Me 


01 


1.06 


11.87 


0.45% 


0.3 - 15 Mj 


209.02 


0.68 


7.62 


3.4% 


25 - 150 Me 


01 


1.05 


11.76 


0.74% 


0.3 - 15 Mj 


877.01 


0.23 


2.63 


8.6% 


6 -AO Me 


02 


0.21 


2.34 


11% 


5-25 Me 


896.02 


0.28 


3.14 


5.9% 


9-30 Me 


01 


0.38 


4.26 


4.9% 


10 - 40 Me 



Note. — The stated fractional uncertainty in the planet radius 
comes from the radius ratio of the planet and star. The stars 
themselves have roughly a 20% fractional uncertainty in their radii 
which has not been included here. The uncertainty in the range of 
planetary masses is dominated by theoretical uncertainties. 
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191.02. 

Next we have KOI 877.01 and 02 have radii of 
2.6Re and 2.3Re, respectively. These planets are 
near the transition to the ice giants Uranus and 
Neptune, but may be volatile-rich sub-Neptunes 



or super-Earths like GJ1214b (Charbonneau et al. 



2009). The estimated mass range is 6 — 40Mb for 
KOI 877.01 and 5 - 2bME for KOI 877.02. Since 
the high-mass, high-density limits are difficult to 
explain by existing planet formation scenarios we 
use estimates of 15 and lOM^;, respectively. 

Three others, KOI 152.02/03 and KOI 896.02, 
have sizes similar to Neptune, so we assign them 
15Mb, noting the large possible range of 9 — 30Mb 
and the anticipated transition to planets possess- 
ing larger H/He envelopes. This transition is 
the reason why KOI 896.01 is estimated at only 
Mp = 20Me despite being much larger than KOI 
896.02. The remaining 4 objects, if confirmed, 
are likely gas giant planets. KOI 152.01 {6.7Re) 
and 209.02 {7.6Re) have radii smaller than Sat- 
urn {9Re, 95Me), hence we assign them a smaller 
mass of 60Me- They could be more massive 
than Saturn with large cores (e.g., HD149026b) 
or lower mass objects with small cores that are 
dominated by a H/He envelope; the plausible mass 
range is consequently wide. The other two, KOI 
191.01 and 209.01, have Jupiter sizes and we as- 
sign them Jupiter masses noting the cautionary 
tale of HAT-P-12b and CoRoT-3b above, which 
also have Jupiter sizes. 

5. Dynamical interactions 

In these data, we detect no significant TTV sig- 
nal given the short time baseline. However, using 
the measured periods and estimates of the plane- 



tary masses from Section 4.3 , we conduct a Monte 
Carlo study to determine what TTV signals we ex- 
pect from these systems with more data. For this 
study we assume coplanar orbits; since inclination 
affects the TTV signal at second order, these re- 
sults apply to systems with < 0.1 radians, where 
(j) is the mutual inclination of the planets. 

5.1. Monte Carlo outline 

The masses and periods of the planets used 
for the Monte Carlo study are fixed to the val- 
ues shown in Table [5l Since the short-term TTV 
signal (St) scales in a known manner with planet 



mass {St ^ JTipcrt for non-resonant systems and 
5t mtrans/(™trans + Wpert) for rcsouaut Sys- 
tems) it is straightforward to adjust these results 
for other planetary masses. Here we define 6t to 
be the Root Mean Square (RMS) of the residuals 
after subtracting a linear ephemeris. 

Parameters that were adjusted in this Monte 
Carlo study include the eccentricity of both plan- 
ets, the longitudes of pericenter, and the mean 
anomalies at the initial time. The eccentricities 
were chosen from a mixture of an exponential and 
a Rayleigh distribution ( Juric and Tremaine|2008 
Zakamska et al.||2010[ ): 



ecc{x) = aXe + (1 



(1) 



where a — 0.38 gives the relative contributions 
of the two probability density functions, A = 15 
is the width parameter of the exponential distri- 
bution, and (Je — 0.17 is the scale parameter of 
the Rayleigh distribution. The value of ae was 
found by fitting the distribution of eccentricities 
in known multiplanet systems measured from ra- 
dial velocity surveys using only systems with mea- 
sured eccentricities. The value of A accounts for 
the planets found to be in or near circular or- 
bits. The values of the longitudes of pericenter and 
mean anomalies were chosen from a uniform dis- 
tribution rather than using the observed relative 
longitudes of the various planets in the systems. 
Thus, these results are applicable to general sys- 
tems with the measured period ratios — including 
the five systems presented here. For each system, 
a large sample of initial conditions was generated. 
Each realization was integrated for seven years, 
the time baseline of an extended Kepler mission. 
The transit times of the planets in the system were 
tabulated and a linear ephemeris was subtracted. 
Finally, the RMS value of the timing residuals was 
recorded. 

A full-scale investigation of the stability of the 
systems used in the study was not feasible due to 
the computational cost. Nevertheless, a few sim- 
ple criteria were used to eliminate systems that 
are most likely to be unstable. First, if any two 
planets came within two Hill radii of each other 
(twice the sum of the Hill radii of the two planets) 
then the system was rejected. Second, if the semi- 
major axis of any planet changed by more than 
20% from its initial value then the system was re- 
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Table 5 

Orbital and physical properties of the systems used in the Monte Carlo study in order 
OF increasing orbital period for each system (not necessarily in order of candidate 

IDENTIFICATION NUMBER). 



System 


Stellar Mass 


Mass 1 


Period 1 


Mass 2 


Period 2 


Mass 3 


Period 3 




(Mo) 




(Days) 




(Days) 




(Days) 


152 


1.4 


15 Me 


13.5 


15 Me 


27.4 


60 Me 


51.9* 


191 


0.9 


10 Me 


2.42 


1 Mj 


15.4 






209 


1.0 


60 Me 


18.8 


1 Mj 


49.3* 






877 


1.0 


15 Me 


5.96 


10 Me 


12.0 






896 


0.8 


15 Me 


6.31 


20 Me 


16.2 







Note. — An asterix indicates an orbital period estimated from the dm'ation of a single 
transit event. 



jected. Third, if the resulting TTV signal for a 
given planet was larger than the period ratio of 
the planets (greater than unity) times the period 
of the planet then the system was rejected ( Agol 
et al.| ( |2005[ ) showed that a TTV signal of order 
the period of the planet is possible, but that oc- 
curs only in the most favorable configuration of a 
2:1 MMR and a very massive perturbing planet). 
We note that this third criterion eliminated only 
a small fraction (<C 1%) of the systems under con- 
sideration. Table [6] indicates the number of orbits 
of the innermost planet used in the study, the ini- 
tial number of systems used in the study, and the 
final number of systems that satisfied the stability 
critera above. 

5.2. Monte Carlo results 

Here we present some of the results from this 
study and give an estimate for the expected TTV 
signal for the five systems considered. We use as 
an example the KOI 896 system and then show the 
essential outcomes for all five systems. Additional 
information about the results of the Monte Carlo 
results is found in the Appendix. 

The KOI 896 system has two Neptune-size plan- 
ets just outside the 5:2 mean motion resonance. 
Figure [9] shows the distribution of TTV signals 
from the simulation in terms of the signal-to-noise 
ratio (SNR) — that is, the ratio of the size of the 
TTV signal to the mean timing uncertainty ob- 



tained from our pipeline measurements of individ- 
ual transit times. The likely minimum size of this 
TTV signal indicates that additional data should 
show such deviations from a constant period ex- 
cept in a limited set of configurations. For ex- 
ample, the fifth percentile of the distribution still 
gives a TTV signal of a few minutes and has an 
SNR near unity. In addition, the fact that it is not 
precisely situated in an MMR might allow this sys- 
tem to be characterized by the analytic methods 



of Nesvorny and Morbidelli ( 2008 1 rather than the 



full numerical simulations needed for resonant or 
near-resonant systems. 

The width of the resulting TTV distribution 
can be understood as follows. The left-most edge 
(smallest TTV signal) in Figure [9] can be esti- 
mated from a circular case, employing equation 



(31) from Agol et al. (2005), extrapolating (and 
simplifying) from the 2:1 MMR. Here we would 
expect TTV signals for the inner planet of order 

St ~ mpert/m*(Pi/P2)Ptrans ~ 30 scc. The TTV 

signal grows substantially with increased eccen- 
tricity. The largest expected TTV signal should 
be of order the period of the planet in question; 
for KOI 896 this is about 10^ seconds. These two 
bounds on the TTV signal are apparent in the re- 
sults of the Monte Carlo simulation. Aside from 
very few realizations that have TTV signals less 
than 10 seconds, the simulation yields results con- 
sistent with the bounds mentioned above. The 
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Table 6 

Information about the number of systems used in the Monte Carlo study. 



System 


Inner Orbits 


Initial Systems 


Final Systems 


152 


185 


200000 


3477 


191 


1033 


15000 


9069 


209 


133 


15000 


12732 


877 


208 


15000 


14701 


896 


392 


15000 


9566 



median expected signal is roughly 2600 seconds for 
the inner planet and 6700 seconds for the outer; 
the ratio of these two values is near their period 
ratio, also as expected. 



Figure 10 shows the distribution of eccentrici- 
ties that survive the stability criteria. Also shown 
is the initial distribution for the eccentricities — 
the solid curve. While the surviving distribution 
of eccentricities tracks the initial distribution, one 
can see that systems with lower eccentricities are 
somewhat more likely to survive than larger ones. 
Note that this figure does not show the final eccen- 
tricities, rather it shows the initial eccentricities of 
the systems that pass all of the stability criteria. 

Table [t] shows the median, fifth, and 95**^ per- 
centiles for the expected TTV signal and SNR for 
each of the candidate planets. Histograms similar 
to figures |9] and [TO] for each system are found in 
the Appendix. With the possible exception of KOI 
191, each of the planets in these five systems will 
likely have observable transit timing variations by 
the end of an extended Kepler mission. For KOI 
191, even if the TTV signal is small it may yet 
be detectable simply because there will be a large 
number of transits (more than 1000) over the du- 
ration of an extended mission which may compen- 
sate for the low signal-to-noise ratio of the TTV 
signal to the transit time uncertainties. 

The primary reason for the small signal in KOI 
191 is the large ratio of orbital periods, exceeding 
6:1. Thus, the TTV signal is weakened signifi- 
cantly. If the outermost planet in KOI 191 were 
to have an eccentric orbit then it would give a pe- 
riodic TTV signal with a period equal to that of 
the outer planet as described in Section 4 of |Agol| 
et al.| ( |2005[ ) (see also [Borkovits et al.||2003[ ). 

For KOI 209, the expected TTV signal for the 
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Fig. 9.— Distribution of the TTV signal for the 
KOI 896 system. The blue histogram is for the 
interior planet and the orange is for the exterior 
planet. Both histograms overlap for much of the 
domain. 




OA a£ 

Fig. 10. — Eccentricity distributions for the sys- 
tems that survived the stability criteria for the 
KOI 896 system. The blue histogram is for the 
interior planet while the orange is for the exterior 
planet. 
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Table 7 

Mean timing precision (ctt), median and quantile TTV signals expected for each planet in 

THE FIVE SYSTEMS. SiGNAL-TO-NOISE RATIOS ARE GIVEN IN PARENTHESES. 



Candidate 




Median TTV 


5% TTV 


95% TTV 




(sec) 


sec/(SNR) 


sec/(SNR) 


sec/(SNR) 


152.03 


790 


2590 (3.3) 


223 (0.3) 


36100 (45) 


02 


700 


39500 (56) 


4090 (5.8) 


387000 (550) 


01 


200 


18900 (97) 


1830 (9.4) 


285000 (1500) 


191.02 


1700 


29 (0.02) 


6 (0.003) 


5850 (3.4) 


01 


70 


20 (0.3) 


9 (0.1) 


3850 (54) 


209.02 


330 


2230 (6.8) 


488 (1.5) 


79600 (240) 


01 


151 


2300 (15) 


234 (1.5) 


73100 (480) 


877.01 


850 


16500 (19) 


2620 (3.1) 


41900 (49) 


02 


760 


79100 (100) 


12500 (16) 


200000 (260) 


896.02 


680 


2630 (3.9) 


103 (0.2) 


80900 (120) 


01 


370 


6700 (18) 


251 (0.7) 


199000 (540) 



inner planet shows an abrupt cutoff and is ex- 
pected to be larger than a few hundred seconds. 
This is because KOI 209.02 has the longest period 
of ah of the inner planets. Given the time basehne 
of the extended mission, planets with periods of 
a few tens of days will likely prove to be among 
the most interesting for TTV studies as they si- 
multaneously have longer periods (the TTV signal 
is linear in the period) and will have a sufficient 
number of transits for a complete analysis. 

The proximity of KOI 877 to the 2:1 MMR in- 
dicates that this system is likely to have very large 
variations. However, a steep drop in the expected 
signal occurs when the orbits are nearly circular. 
Figure [TT] shows an expanded view of the TTV sig- 
nal for the inner planet in KOI 877 as a function 
of the inner and outer planet eccentricities. From 
this figure one can see that, while the zero eccen- 
tricity case exhibits a relatively small TTV signal, 
eccentricities much larger than 0.01 cause the sig- 
nal to increase beyond an SNR of unity near 10^ 
seconds (~ 15 minutes). Should the TTV signal be 
this size or smaller, it should stringently constrain 
the eccentricities of both planets in the absence of 
any other data. We note that all of these results 
for the expected TTV signal have significant de- 
pendence on the eccentricities of the planets. One 
consequence of this fact is that, if a large fraction 
of these or other multiple systems do not show a 



TTV signal, then low eccentricity orbits are much 
more common in multi-planet systems than in sin- 
gle planet systems. 

The three-planet system of KOI 152 portends 
the exciting and challenging studies of systems 
where there are more than two planets and where 
multiple planets transit the star. This system is 
particularly interesting given the relatively close 
proximity to the 4:2:1 multibody resonance. How- 
ever, it is unlikely that this system occupies this 
resonance given the estimated orbital periods of 
the planets — one being estimated from a single 
transit. For KOI 152, the middle planet is likely to 
exhibit the largest TTV signal — being just outside 
the 2:1 MMR with an interior planet and perhaps 
just interior to the 2:1 MMR of the exterior planet. 

5.3. Challenges for systems vifith more 
than two planets 

One challenge that three-planet systems, such 
as KOI 152, pose is the confusion that can arise 
from multiple, competing perturbers in the TTV 
signal for a particular planet. We present three 
broad scenarios for consideration in future stud- 
ies, although other regimes may exist: 1) nonreso- 
nant/nonresonant where there is no mean motion 
commensurability between any pair of planets, 2) 
resonant /nonresonant where one pair of planets 
has a mean motion commensurability while the 
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Fig. 11.— Contour plot of the TTV signal for the 
inner planet in KOI 877 as a function of the in- 
ner and outer planet eccentricities. Note that this 
is an expanded view of the lower-eccentricity sys- 
tems, the eccentricity distributions of both plan- 
ets extend beyond 0.6. The contours correspond 
to signal-to-noise ratios of 1 and 10 (times of 850 
and 8500 seconds). 



other does not, and 3) resonant/resonant where 
any pair of planets lies near a mean motion com- 
mensurability. 

For the first scenario, the TTV signal due to 
one perturber should be largely independent of 
the TTV signal due to the second perturber. The 
effect from both perturbers will be of order the 
perturber to stellar mass ratio, and therefore may 
be comparable. But their contributions will con- 
tribute linearly to the overall signal and the peri- 
odicities in the TTV signal due to one perturber 
will be independent of the periodicities induced 
by the other. In other words, a Fourier transform 
of the TTV signal would likely show two sets of 
independent peaks (see Steffen|[2006 1 that can be 
distinguished provided the data have a sufficient 
time baseline (iobs ^ 1/^/ where A/ is the typi- 
cal difference in frequency of the most prominent 
Fourier components between the TTV signals of 
each planet). 

One's ability to identify the orbital elements of 
the planets in the second scenario, with a reso- 
nant pair and a nonresonant companion, will de- 
pend upon which planet is transiting. If a planet 
in the resonant pair is transiting, then the system 
is likely to be invertible for two reasons. First, 
the resonant signal will be significantly larger (by 
a factor of the stellar mass to the total planetary 
mass of the resonating system) and will therefore 
be more readily identified due to its enhanced sig- 
nal to noise ratio. Second, similar to the first sce- 
nario, the TTV contribution from the nonresonant 
perturber will combine linearly with those of the 
resonant perturber and sufficient data will distin- 
guish their contributions. 

Characterizing a system where the nonresonant 
planet is transiting may be much more challeng- 
ing. Here, the transit times will vary on timescales 
of the libration period of the perturbing, resonat- 
ing pair of planets. In general, a signal with a sim- 
ilar period and amplitude may also be generated 
with a single planet whose orbital period is equal 
to a multiple of the libration period of the res- 
onating pair. In addition, resonating systems may 
exhibit secular evolution on timescales of only a 
few years, which could also affect the nonresonant 
planet in a measurable way. 

The third scenario, where the system has mul- 
tiple pairs of planets near MMR, may present seri- 
ous challenges if only one planet transits. In favor- 
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able configurations, sufficient data will allow one 
to identify the dynamical interactions among var- 
ious planets. This may be more feasible when the 
system architecture, while resonant, is also hier- 
archical. Consider, for example, a 1:2:8 hierarchy 
where the 4:1 resonance between the outer two 
planets is likely to be much stronger than the as- 
sociated 8:1 resonance between the outermost and 
innermost planets. On the other hand, compact 
and strongly interacting systems such as 1:2:4 or 
1:2:3 may produce TTV signals that mimic single 
perturbing systems that are in a different reso- 
nance altogether (e.g., a 3:2 or 4:3 MMR). 

All of the challenges in characterizing multi- 
body systems from the TTV signal that are listed 
above are lessened significantly when multiple 
bodies transit the star. Most importantly, in a 
multiply transiting system the periods of more 
than one planet are known — eliminating confu- 
sion regarding, for example, which resonance some 
planets may occupy. In addition, if there is a non- 
transiting perturbing planet in a system with two 
or more transiting planets then its effects will be 
present, though perhaps very small, in the TTV 
signal of each of the transiting planets. Thus, cor- 
rectly identifying non-transiting perturbing plan- 
ets will likely be easier in a multiply transiting 
system, such as those presented here, than in a 
singly-transiting system. 

We illustrate this point with a realization of the 
KOI 152 system, given by parameters in Table |8] 
starting with circular orbits. Figure 12 shows the 



transit timing variations for all three planets and 
for each pair of two planets separately. If the in- 
ner planet failed to transit, it may be difficult to 
detect from the transit times of the other planets 
(compare panels a and b). However, if the outer 
planet failed to transit, its presence would be be- 
trayed because the two inner, transiting planets 
would not oscillate with the same frequency, in 
antiphase from each other, as they would if they 
were alone (panel c) . Finally, if the middle planet 
failed to transit, the large TTVs of both the inner 
and outer planet would be too large to be due to 
each other, and the economical hypothesis of an 
intervening planet could explain both their TTV 
patterns. 




200 400 6D0 BOO 1000 1200 
t- 24549 00 

Fig. 12. — Simulated transit times, relative to a 
linear ephemeris for a system based upon KOI 
152. (a) The full, 3-planet system with parame- 
ters given by Tablejsj (b) The same system, except 
the inner planet is absent. The two outer planets 
interact very similarly to the case with all three 
planets, (c) Now the exterior planet is absent. 
The inner planet interacts with the middle planet 
much as in the full system, but the middle planet, 
now with the dominant interaction of the outer 
planet absent, oscillates in antiphase with the in- 
ner planet, (d) Now the middle planet is absent. 
The inner and outer planets are too widely sepa- 
rated for a significant TTV signal to be present. 
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Table 8 

Integration parameters for the KOI 152 representative system shown in Figure [T2l 



Candidate 


P (Days) 


A (deg) 


e 


771 (Earth) 


03 


13.48 


-12.551885 


0.00 


15 


02 


27.40 


71.211695 


0.00 


15 


01 


51.94 


-74.197154 


0.00 


60 



Note. — The initial epoch is t-2454900 = 120 and the 
assumed stehar mass is 1.4 Mq 



6. Discussion 

We presented five Kepler targets, each of 
which has muhiplc transiting exoplanet candi- 
dates. These candidates have not been fully vet- 
ted and therefore none is a confirmed exoplanet. 
Yet, each of these systems have passed several im- 
portant validation steps that are used to eliminate 
false positive scenarios. It is difficult to construct 
viable astrophysical solutions that are consistent 
with all of the data. Of particular importance 
is the fact that many common false positives for 
singly transiting systems are not viable false pos- 
itives for multiply transiting systems. For ex- 
ample, dynamical stability precludes triple star 
systems with orbital period ratios of order unity. 
Another possibility is a foreground star with two 
background eclipsing binaries, which requires two 
background eclipsing binaries to be blended within 
the same point-spread- function (PSF). Finally, it 
is unlikely that two systems that have periodic 
planetary transit features will have period ratios 
very near two-to-one. 

Additional data, such as high precision radial 
velocity measurements, may provide sufficient ev- 
idence to confirm these as planets or refute them 
as instances of a false positive signal. At this time 
these systems are not a sufficiently high priority 
for the Kepler team to conduct such observations. 
If these systems are planets and they follow the ob- 
served mass-radius relationship of known planets, 
the Kepler mission is likely to find TTV signatures 
due to their mutual interactions with the possible 
exception of KOI 191. Over the course of the mis- 
sion, a detailed TTV analysis of these systems can 



constrain their libration amplitudes, masses, and 
eccentricities — confirming them as planets with- 
out the need for many high-precision RV measure- 
ments (as in Kepler-9 ( Holman et aL||2010 )). Ad- 
ditionally, TTV measurements from the transits 
of multiple planets can provide better estimates of 
the stellar properties (e.g., density, limb darkening 
model) than systems with only a single transiting 
planet. 

The possible observation of Transit Duration 
Variations (TDVs) in these systems may identify 
orbital precession, moons, or Trojan companions 



(Ford and Holman 2007 Kipping 2009 Kipping 



et al. 2009). For some multiple transiting systems. 



it should be possible to constrain the relative or- 
bital inclination based on TTVs, TDVs, and the 
constraint of orbital stability. Once a sizable and 
minimally-biased population of multiple transiting 
systems has been identified (rather than this small 
and select sample), it will also be possible to char- 
acterize the frequency of multiple planet systems, 
the distribution of mutual orbital inclinations, and 
identify the orbital architectures of planetary sys- 
tems. Collectively, this information will provide 
considerable insight into the formation, migration, 
and dynamical evolution of planetary systems. 
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A. Appendix 



In this appendix we present the balance of the results from the Monte Carlo study of the TTV effect for 
these systems. In particular, we start with KOI 877, then 209, 191, and finally the three-candidate system 
of 152. 



A.l. KOI 877 



The KOI 877 (Figure 13) system lies very near, or possibly within, the 2:1 MMR; and therefore will 
likely have a large TTV signal. The Monte Carlo study shows this as the median TTV signal from this 
system is several hours. Interestingly, the systems that survive the stability criteria have a virtually identical 
eccentricity distribution to the initial distribution. Given the proximity to MMR, it is likely that an inde- 
pendent, in-depth study of the long-term dynamics would reject more of the high eccentricity systems and 
the expected TTV signal would decline somewhat. Nevertheless, a system with nearly equal masses near a 
2:1 MMR is an ideal scenario to find a TTV signal — even at zero eccentricity. 



A.2. KOI 209 

For KOI 209, the smaller inner planet is likely to show the largest TTV signal. Interestingly, the surviving 
eccentricity distribution of the more massive outer planet is very skewed toward smaller values while the 
eccentricity of the inner planet more closely matches the initial distribution. The surviving eccentricity 



distributions of the KOI 209 system (Figure 14) strongly favor a low-eccentricity outer planet. Consequently, 
the expected TTV signal for the outer planet has a tail toward lower values. This is typical of systems where 
the planets are not near a mean motion resonance, but interact on secular timescales to exchange significant 
angular momentum. The larger planet mass and semimajor axis means that a similar eccentricity results 
in a much larger angular momentum deficit for the outer planet. The maximum eccentricity of the inner 
planet over a secular timescale can be more sensitive to the initial eccentricity of the outer planet than its 
own initial eccentricity. 



A.3. KOI 191 



The hierarchical structure of the KOI 191 system (Figure 15 ) means that a TTV signal is likely to be 
significant only if there is significant eccentricity in the orbit of either planet. The sharp peak in the TTV 
distribution of the outer planet near 10 seconds agrees with the expectation for an outer transiting planet 
in a system where the planet interactions are negligible (section 3 of Agol et al. (2005 )). In this scenario the 



outer planet probes the astrometric deviations of the star due to the inner planet. A rough estimate of the 
size of this effect can be found by taking the characteristic displacement of the star from the barycenter and 
divide by the characteristic velocity of the outer planet (aimi/m»)/(a2/i-2) — lOsec. It may be possible to 
observe variations in the duration of the transit of the outer planet due to the movement of the star within 
the inner binary while the outer planet is transiting in other hierarchical systems, but perhaps not with the 
KOI 191 system. 



A.4. KOI 152 

For this system, all three planets should have a sizeable TTV signal due to the proximity of the 2:1 
MMR between the inner and outer pair. As mentioned above, the middle planet will likely show the largest 
TTV signal. One would expect the TTV signal for the outer planet to be roughly half that of the middle 
planet — it has a factor of two larger period, but the perturbing, middle planet is only one quarter the mass. 
The eccentricity distribution of the outer two planets favor smaller eccentricities. 
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SuniiriiigamihmlT TTVSNR 

Fig. 13. — Left: Distributions of eccentricities that survive the stabiUty criteria for the KOI 877 system. 
Right: Distributions of the TTV signal for the KOI 877 system. The blue histogram is for the interior planet 
while the orange is for the exterior planet. 




SuniiriiigamihmlT TTVSNR 

Fig. 14. — Left: Distributions of eccentricities that survive the stability criteria for the KOI 209 system. 
Right: Distributions of the TTV signal for the KOI 209 system. The blue histogram is for the interior planet 
while the orange is for the exterior planet. 
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SuniiriiigamihmlT TTVSNR 

Fig. 15. — Left: Distributions of eccentricities that survive the stability criteria for the KOI 191 system. 
Right: Distributions of the TTV signal for the KOI 191 system. The blue histogram is for the interior planet 
while the orange is for the exterior planet. 
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Fig. 16. — Left: Distributions of eccentricities that survive the stabihty criteria for the KOI 152 system. 
Right: Distributions of the TTV signal for the KOI 152 system. The blue histogram is for the interior planet, 
orange is for the middle planet, and light gray is for the exterior planet. 
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